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Dwarf spiral galaxies - and in particular the prototypical DDO 154 - are known to be completely 
dominated by an unseen component. The putative neutralinos - so far the favored explanation for 
the astronomical dark matter - fail to reproduce the well measured rotation curves of those systems 
because these species tend to form a central cusp whose presence is not supported by observation. 
We have considered here a self-coupled charged scalar field as an alternative to neutralinos and 
investigated whether a Bose condensate of that field could account for the dark matter inside DDO 
154 and more generally inside dwarf spirals. The size of the condensate turns out to be precisely 
determined by the scalar mass m and self-coupling A of the field. We find actually that for m 4 /A ~ 
50 - 75 eV 4 , the agreement with the measurements of the circular speed of DDO 154 is impressive 
whereas it lessens for larger systems. The cosmological behavior of the field is also found to be 
consistent - yet marginally - with the limits set by BBN on the effective number of neutrino families. 
We conclude that classical configurations of a scalar and self-coupled field provide a possible solution 
to the astronomical dark matter problem and we suggest further directions of research. 



I. INTRODUCTION 

After many years of global consensus on the fact that dark matter consists in Weakly Interacting Massive Particles 
(WIMPs) - like, for instance, the lightest neutralino in the Minimal Supersymmetric Standard Model (MSSM) - there 
is still no strong evidence in favor of WIMPS, neither from bolometer experiments designed for direct detection, nor 
from the observation of cosmic rays, a fraction of which could consist in WIMPs annihilation products. This absence 
of experimental constraints on dark matter from the particle physics side leaves the door wide open for alternative 
descriptions of the hidden mass of the Universe. 

Moreover, in the past three years, there has been a lot of controversy concerning the small-scale inhomogeneities 
of the WIMPs density. Indeed, many recent N-body simulations of structure formation in the Universe suggested 
that any dark matter component modelized as a gas of free particles - such as WIMPs - tends to cluster excessively 
on scales of order 1 kpc and smaller. This would result in cuspy density profiles at galactic centers, while most 
rotation curves indicate a smooth core density [1]. Many galaxies even seem to be dominated by baryons near their 
center, with a significant dark matter fraction only at large radii. This clearly contradicts the results from current 
N-body simulations, in which the dark matter density is strongly enhanced at the center of the halo with respect to 
its outskirts. 

This argument was attacked by Weinberg and Katz [2] , who stressed the importance of including the baryon component 
in N-body simulations. Indeed, the baryon dissipative effects could be responsible for a smoothing of the central dark 
matter cusp in the early Universe. This possible solution to the dark matter crisis was discarded later by Sellwood 
[3], who found opposite results in his simulation. 

Apart from the central cusp problem, N-body simulations raised some secondary issues [1]. First, a clumpy halo 
could generate some tidal effects that could break the spatial coherence of the disk. Second, the predicted number of 
satellite galaxies around each galactic halo is far beyond what we see around the Milky Way. Third, the dynamical 
friction between dark matter particles and baryons should freeze out the spinning motion of baryonic bars in barred 
galaxies. All these arguments are still unclear, because they seem to depend on the resolution under which simulations 
are carried [4,5], and also because of our ignorance of what could be the light -to-mass ratio inside small dark matter 
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clumps. In addition, the predicted number of satellite galaxies doesn't seem to be in contradiction with constraints 
from microlcnsing [6]. 

Should these various problems be confirmed or not, it sounds reasonable to explore alternatives to the WIMPs model, 
or more generally speaking, to any description based on a gas of free particles. This can be done in various ways: for 
instance, one can introduce some deviations from a perfect thermal phase-space distribution [7] , or add a self-coupling 
between dark matter particles [8] . A more radical possibility is to drop the assumption that dark matter is governed 
by the laws of statistical thermodynamics. This would be the case if dark matter consisted in a classical scalar field, 
coherent on very large scales, and governed by the Klein-Gordon equation of motion. 

This framework should be clearly distinguished from other models of bosonic dark matter, like those based on heavy 
bosons - for instance, sneutrinos - or axions. In the first case, the Compton wavelength h/(mc) of an individual 
particle is much smaller than the typical interparticle distance while in the second case, for axion masses of order 
10~ 6 eV, it is still much smaller than the typical size of a galaxy. So, in these examples, the bosons can be described 
on astronomical scales like a gas of free particles in statistical equilibrium. It follows that the halo structure cannot 
be distinguished from that of standard WIMPs. 

A coherent scalar field configuration governed by the Klein-Gordon and Einstein equations is nothing but a self- 
gravitating Bosc condensate. Such condensates span over scales comparable to the Dc Broglie wavelength, L = h/p. 
In the case of free bosons - i.e., with a quadratic scalar potential - the momentum p is of order mv csc where v esc 
is the escape velocity from the system. Typical examples are boson stars [9-12], for which the characteristic orders 
of magnitude discussed in the literature are, for instance, m ~ 10 GeV and w osc ~ c, leading to a radius as tiny as 
L ~ 10~ 14 cm. Even for axions, which have a much smaller mass and an escape velocity given by the motion of stars 
in a galaxy - v csc ~ 100 km/s - the Dc Broglie wavelength is only of order L = 100 km, so that on galactic scales, 
the medium can be treated as a gas. In order to obtain a galactic halo described by the Klein-Gordon equation, one 
should consider masses of order m = h/(Lv csc ) where v csc ~ 100 km/s and L ~ 10 kpc. This yields m ~ 10 -23 eV. 
Such an ultra-light scalar field was called "fuzzy dark matter" by Hu [13], who discussed its overall cosmological 
behavior. In some previous works, we focused on a variant of this model in which the ultra-light scalar field is 
complex - then, the conserved number associated with the U(l) global symmetry helps in stabilizing the condensate 
against fragmentation [10,12]. In [14] - thereafter Paper I - we compared the rotation curves predicted by this model 
with some data from spiral galaxies. In [15] - thereafter Paper II - we simulated the cosmological evolution of the 
homogeneous background of such a field. The model seems to be quite successful in explaining the rotation curves, 
but it has two caveats. First, such a low mass is very difficult to implement in realistic particle physics models. 
The second problem is related to the fact that because of the U(l) symmetry, the field carries a conserved quantum 
number. As explained in Paper II, the value derived from cosmological considerations for the density of this quantum 
number does not seem consistent with that inferred from astrophysical arguments. 

These caveats motivate the introduction of a quartic self-coupling term in the scalar potential. In that case, it is 
already known from boson stars that for the same value of the mass, the self-coupling constrains the field to condensate 
on much larger scales [16,12] - the size of the self-gravitating configurations is still given by L = h/p, but in presence 
of a self-coupling, the momentum cannot be identified with mv csc . Then, without recurring to ultra-light masses, we 
may still describe the galactic halos with a Bose condensate. A massive scalar field with quartic - or close to quartic 

- self-coupling was proposed as a possible dark matter candidate by Peebles, who called it "fluid dark matter" [17]. 
In this paper, we will study a variant of fluid dark matter in which the quartically self-coupled massive scalar field is 
complex, still for stability reasons 1 . Some cosmological properties of such a field were already discussed in [19]. 

We will focus mainly on galaxy rotation curves, assuming that the dark matter halos are the self-gravitating, aspherical 
and stable equilibrium configurations of our scalar field in the presence of a baryonic matter distribution - stellar disk, 
HI gas, etc... We will present here the first solution of this problem. However, we should stress that some different 
models in which the rotation curves are also seeded by a coherent scalar field were studied previously by Schunk [20] 

- with a vanishing scalar potential, Goodman [21] - with a repulsive self-interaction, Matos et al. [22], Nucamendi et 
al. [23], Wetterich [24], Urena-Lopez and Liddle [25]. In some of these papers, and also in many other recent proposals 

- see for instance [26] the main goal is to try to solve simultaneously the dark energy and dark matter problems, 
assuming that a quintessence field can cluster on galactic scales. This raises some subtle issues, like the existence of 
a scalc-dcpcndcnt equation of state. At the present stage, we do not have such an ambition, and we will focus only 
on the dark matter problem. 

In section II, we write the Einstein and the Klein-Gordon equations which govern the scalar field and the gravitational 
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potential distributions in the presence of a given baryonic matter density. We will see that these equations can be 
combined into a single non linear Poisson equation. The solutions are technically difficult to find, first, due to the 
non-linearity, and second, because some boundary conditions are given at the center, some others at infinity. So, it is 
not possible to follow a lattice approach, in which one would start from a particular point and integrate numerically 
grid point by grid point. However we present in section III a recursive method which allows to find all the exact 
solutions after a few iterations. In section IV, we compare the galaxy rotation curves obtained in this way with some 
observational data. We lay a particular emphasis on the dwarf spiral galaxy DDO 154, for which the rotation curve is 
among the most difficult to explain with usual dark matter profiles. We will see that a mass-over-self-coupling ratio 
to 4 / A ~ 50 (cV) 4 provides a very good fit to the DDO 154 rotation curve - however at the expense of poor fits to 
the largest spiral galaxies. Because the scalar field condensates inside the gravitational potential wells of baryons and 
strengthens them, the question of its effects on the inner dynamics of the solar system naturally arises. We derive in 
section V the modification of the solar attraction in the presence of the self-interacting scalar field under scrutiny and 
show that an anomalous acceleration appears that is constant and that points towards the Sun. We investigate the 
limit set on our model by the Pioneer radio data. In Paper II, we studied the cosmological behavior of a homogeneous 
scalar field that was assumed to play the role of dark matter at least from the time of matter-radiation equality until 
today. This analysis is updated in section VI where we specifically assume to 4 / A ~ 50 (eV) 4 . Such a large value point 
towards a large total density of the Universe during radiation domination which is at the edge of the current bounds 
set in particular by BBN on cosmological parameters. The last section is devoted to a discussion of the strong and 
weak aspects of our alternative dark matter model. We finally suggest some further directions of investigation beyond 
the simple but restrictive framework of isolated bosonic configurations. 



II. GRAVITATIONAL BEHAVIOR 



The complex scalar field cf> under scrutiny in this article is associated to the Lagrangian density 

£ = g^d^d^ -V(4>) , (l) 

where the U(l) invariant potential V includes both quadratic and quartic contributions 

V{4>) = toV^ + A {<^} 2 . (2) 

The gravitational behavior of the system follows the standard GR equations whilst the field <j> satisfies the Klein- 
Gordon equation 

d»{^—gg^d^} + f£ = , (3) 



where g^ u denotes the metric. We would like to investigate to which extent the scalar field 4> m ay account for the 
dark matter inside galaxies. The problem simplifies insofar as the gravitational fields at stake are weak and static. 
In this quasi-Newtonian limit of general relativity, deviations from the Minkowski metric rj^ = diag{\, —1, — 1, — 1} 
are accounted for by the perturbation tensor - from now on, we use the convention c = 1. The Newtonian 
gravitational potential <f> = h on /2 is actually a small quantity of order v 2 sc ~ 10~ 7 — 10 -6 , where u csc denotes the 
escape velocity. Our analysis is based on an expansion up to first order in <t>. The baryonic content of galaxies is 
described through the energy-momentum tensor 

T"" = {p h +P h )U»U v - P h9fW , (4) 

where [/ M = {1,?7}. Baryons behave as dust with non-relativistic velocities. Actually, because galaxies are virialized 
systems - hence the assumption of static gravitational fields - the spatial velocity v is a small quantity of order 
Vcsc ~ V 7 ^. The kinetic pressure-to-mass density ratio is even more negligible since Pb/pb ~ v 2 ~ v 2 BC ~$<1. We 
are interested in classical configurations where the field (/) is in a coherent state such as 

(j ( 

4>{x, t} = — ^ exp (— iujt) . (5) 

Indeed, one can prove that all stable spherically symmetric configurations can be parameterized in that way [10]. The 
time-derivative <9o</> equals —iuxf), whereas the space-derivative di<f> is of order <p/L where L is the physical length of 
the configuration. That length - which is related to the parameters to and A of the potential - is required to be ~ 1 
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- 100 kpc to account for the galactic dark matter. On the other hand, we shall see later that u is very close to the 
mass to, which is numerically found to be in the ballpark of a fraction of eV. We readily infer a ratio 

I* ' 6 .4X10- mm . (6) 

So, the space derivative di(j> of the scalar field can be safely neglected throughout the analysis. 

In its weak-field limit, general relativity becomes a gauge theory. By conveniently choosing the gauge of harmonic 
coordinates in which the metric perturbation h^ lL , satisfies the condition 

d»h» a = \d a {h^} , (7) 

the GR equations simplify into 

□ V = -WttGS^ . (8) 
The effective source S^ v is related to the energy-momentum tensor T^ v through 

S/tv = - -gp V T\ . (9) 

In the propagation eq. (8), the source is computed in flat space while the metric perturbation is of order $. 
If the dark matter inside galaxies is understood as some classical configuration of the field <£>, should take into 
account both baryonic population - stars and gas - and scalar condensate. In the Newtonian limit where gravito- 
magnetic effects are disregarded, the non-relativistic velocities of baryons can be neglected. The only non-vanishing 
components of the baryonic source tensor are 

S h o = Y and S£ = . (10) 

Assuming that eq. (5) describes the scalar field configuration and disregarding the space-derivatives di<p leads to the 
source components 



where the potential is 



S& = wV - V whilst S t J = -rnjV , (11) 



V{a) = ^a 2 + \a± . (12) 



The well-known solution of the Lienard and Wiechert retarded potentials satisfies the propagation eq. (8). We readily 
conclude that the metric does not contain any space-time component hoi and may be expressed at this stage as 

dr 2 = (1 + 2<$>)dt 2 - (1 - 2^>)5 lJ dx l dx j , (13) 

where the static potentials $ and W are given by integrals over the source distribution V of the baryonic and scalar 
mass densities 

m = -G [ J2-r { Ph (y) + p'4y)\ , (14) 
jt> f — y\ 1 J 

and 

*(f) = —G f -S^ {p b (y) + 4(y)} ■ (15) 
Jv \x - y\ t v > 

The densities p'^ and p'L are respectively defined by 

P ; = 2,V-mV-^ 4 , (16) 

and 
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p£ = mV 2 + ^. (17) 

The potentials $ and \& are different a priori. A careful inspection of the Klein Gordon equation will eventually show 
that they are actually equal. The latter may be written as 

(1 - 2$) 4> - (1 + $ - 3*)- 1 d t {(1 + $ - dicj)} + J^l = . (18) 

The space-dependent term is some 53 orders of magnitude smaller than its time -dependent counterpart and we can 
safely disregard it so that relation (18) simplifies into 

Act 2 = (1 - 2$)w 2 - to 2 , (19) 

where the configuration (5) has been assumed. The scalar field is in a classical state that may be pictured as a Bose 
condensate on the boundaries of which the gravitational potential is 

--K'-£) ■ (20) 

Because the potential $o is a small quantity, the pulsation u> is very close to the mass m. The scalar field essentially 
vanishes outside the condensate whereas its inner value is directly related to the gravitational potential <& through 

Act 2 = 2 ($ - $) w 2 ~ 2 ($ - $) m 2 . (21) 

This relation has important consequences. To commence, the densities p'^ = m 2 a 2 {1 + 3$o — $} and p"^ = 
m 2 a 2 {1 + 3>o — — 2'f} become both equal to p$ ~ m 2 a 2 at lowest order in the potentials. Then $ = \t and 
the metric simplifies. It is straightforward to show that it readily satisfies the gauge condition (7). The scalar field 
density may be expressed as a difference between the gravitational potentials inside and on the boundary of the scalar 
field condensate 

2 to 4 

P4> = ( $ o-$) H(*o-*) , (22) 

where Tt(x) = 1 for x > and Ti(x) = elsewhere. This leads to the Poisson equation 

4 

m 

A$ = 4irG Ph + 87rG— ($ - $) W($ - *) • (23) 

A 

Inside the condensate, gravity turns out to be effectively modified by the presence of the scalar density p<f, whereas 
the conventional Poisson equation is recovered outside. Defining the Planck mass through Mp = l/y/G, we derive a 
typical scale of 

2 A Mp A 

for the scalar field configurations in which we are interested. The dimcnsionless constant A has been introduced by 
[16] in their analysis of self-interacting boson stars. The scale L is related to the mass m and the quartic coupling A 
through 

so that values of the mass in the ballpark of the eV may well be compatible with a size L of order a few kiloparsecs. 
As already noticed by [16], the space-dependent term in the Klein-Gordon eq. (18) is actually suppressed by a factor 
of A which, in our case, reaches values as large as ~ 10 53 . The key feature of the scalar field configurations at stake 
is the existence of a unique scale L that depends only on the parameters m and A of the potential V. 

A pure scalar field configuration may also be seen as a mere fluid with mass density p^. The pressure may be 
derived from the space-space component = — r)ij P$ of its energy-momentum tensor. This leads to 

= C = g 0Q <pU - ^ - X ( $ o-1>) 2 (26) 
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inside the condensate where $ < <£>o- The corresponding equation of state boils down to 

p* = Aw 2 ■ ( 27 ) 

4 771 

and features the generic polytropic form P = K p r where T = 1 + 1/n and isT = A/4m 4 are constants. Spherical 
symmetric solutions of configurations in hydrostatic equilibrium are searched of the form pj pc = & n (z) and P/ Pq = 
0' l+1 (z) whgj-Q z = r/L is the dimensionless radius. The typical scale of the polytropc depends on the central density 
pc and pressure Pq through 

^ = ^£ . (28) 
whereas the generic function satisfies the Lane-Emdcn equation 

with the initial conditions 0(0) = 1 and 0'(O) = 0. In the scalar field case, the polytropic index is n = 1 and the 
solution 

= e(z) = ™£ (so) 

P<t>c z 

readily obtains. It describes a spherical symmetric configuration where the scalar field alone is bound by its own 
gravity. The radius of the pure scalar field condensate is then R = it L where the scale L has already been derived in 
relations (24) and (25): 

k i 1/2 r a i 1/2 i 



R = ttL = tx I } =%i } — (31) 

\2ttGJ \8ttGJ m 2 V ; 

The effect of an asphcrical distribution of baryons on the scalar field condensate will be examined in the next section. 

III. RESOLUTION METHOD 

We would like to compute the gravitational potential $ associated with any density of baryons pb in the galaxy, in 
the presence of a scalar field condensate. So, we need to solve eq. (23). The Heaviside function renders this equation 
strongly non-linear: different solutions have different surfaces where 

$(£ ) = $ , (32) 

so the sum of two solutions is not a solution. Nevertheless, it is possible to solve the equation with a recursive method. 
The idea is to start from an approximate solution <fr(°), and to find from the iterations 

4 

A$(™ +1 ) = 47rG/3 fc + 8ttG^- ($ - 7i($ -$ (n) ) • (33) 

If, for a judicious choice of the $( ra )'s converge towards a limit <I>( 00 ), then the latter will be an exact solution of 
(23). 

We will always work in the approximation in which the baryonic density Pb{x) is axially symmetric, continuous and 
vanishing at infinity. So, the induced gravitational potential should be 

axially symmetric, 

continuous and twice derivable, (34) 
vanishing at infinity. 

We introduce a spherical coordinate system (r, 9, <p) where (9 = 0) defines the symmetry axis. So, there will be no 
(^-dependence in the solutions, and a good way to find them is to perform a Legendre transformation. Let's first 
illustrate this for the general Poisson equation: 



G 



A$(r, 6) = S(r, 6) , (35) 

where S and <I> possess the properties of symmetry and continuity listed previously. If one decomposes the potential 
and the source term S into Legendre polynomials: 



$(r,6) =^P ; (cos0)$ ; (r), (36) 

1=0 

S(r,6)=J2 p i( cose ) s t(r), (37) 



then the S^'s are found from 



21 + 1 



-l 



Si(r) = ^—J S(r,9)P l (cos9)d(cos9) , (38) 



while the <&z's are the solutions of the linear set of equations 



The boundary conditions are given by the properties (34) for all Us, 

-f-$z(0) = and lim $;(r)=0. (40) 

dr r— >+oo 

So, in order to find the solution of eq. (39), one can first compute some Green functions Gi that are continuous, null 
at infinity, with zero derivative at the center, and verifying 



1 d ( 2d 1(1 + 1), 



where 5 is the Dirac function. The unique answer is 

21 + 1 I W v ' Vu 



G <(^ «) = -*TTT ( & H ( u r ) + ( L ) H ( r - u ) \ ■ ( 42 ) 



r + oo 

Since $;(r) = / Si(u)Gi(u,r)du, one finally finds 
Jo 



Hr) = _/ S^y +2 ^ - — Jf S^ 1 ^. (43) 

We can still use this Green function technique in our recursive method. Indeed, if shares the properties (34), 
then we can identify the right-hand side of eq. (33) with S(r,9), and find using the method described above. 

Then also shares the properties (34). 

So, at each recursion step, we need to expand the right hand-side of eq. (33) in Legendre coefficients. Note that the 
$| n ^'s are known from the previous iteration, while the number <f>o has to be imposed in some arbitrary way. In fact, 
looking again at eq. (23), it is clear that there should be different solutions $ associated with different values of the 
free parameter $o- Intuitively, this parameter tunes the size of the bosonic halo, since it defines the surface inside of 
which the scalar field plays a role. In the recursion technique, a possible strategy could be to impose $o once and for 
all. Proceeding in that way, we found that the solution did not converge properly. In fact, it is much more efficient to 
choose arbitrarily a point of coordinates (ro, 9q), and to impose step by step that this point remains on the boundary; 
in other words, for each n, we define $o a s & n '(ro, 9$). For a given 9o, the different possible choices of ro generate a 
one-parameter family of solutions. However, only a finite range of ro values lead to a solution, from r m j n = for no 
bosonic halo, to 



~ J = n L (44) 



ItGm 

for a pure scalar field configuration - see eq. (31) or the alternative derivation in Appendix A. The choice of 9q itself 

is irrelevant and we checked that any other choice gives the same family of solutions. 

In summary, for a given baryonic density, the modified non-linear Poisson can be solved by: 
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1. choosing an arbitrary 6q. 

2. choosing a value ro reflecting the size of the scalar held halo. 

3. defining a starting function . 

4. integrating eq. (33) recursively, with the help of equations (38), (43) and (36). 

We show in Appendix B how to define a function <&(°> which is close enough to the real solution in order to ensure 
fast convergence. 

Let us illustrate this technique with a particular example, based on a stellar disk plus a bosonic halo. In the following, 
we will always treat the stellar disk as a thin distribution with exponentially decreasing density. The optical radius 
is defined in such a way that it encompasses 83% of the total stellar mass, and the thickness of the disk is chosen to 
be twenty times smaller than its radius, so that 



Pb(r, 9) oc exp I — 3.2 1 cos#|— — > cxp < — 3.27I sin#|— — 

r op t ) I r opt 



(45) 



where 7 = 20. Let us choose a case where tq and r opt are comparable, so that baryons and bosons both have an 
influence, for instance: 



'opt 



and 



(r o ,0 o ) = (2- 



(46) 



We plot on figure 1 the function 



$(")(r, tt/2) 



(47) 



$(™)(r ,7r/ 2 ) 

for n — 0, 1,5, 10. In this example, the value <f> = $0 is reached at r = 2r opt in the disk plane, and r = 1.85 r opt in 
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FIG. 1. Gravitational potential $' n ' seen in the direction 6 = n/2 after n = (in red), 1 (in green), 5 (in blue) and 10 (in 
purple) recursions, with r op t = J"max/7r, ro = 2r ma x/vr, 80 = n/2. The functions are seen to converge quickly. 

the orthogonal direction. The oblate form of the cquipotcntials is seen on figure 2. 
A good test of the recursive method is to pick up different values of 9q, and to see whether there is always a value 
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0.5 1 1.5 2 

x/r 

FIG. 2. Equipotential lines generated by a disk-shaped density and a scalar field halo (in red), compared to spherical 
equipotentials (in blue). This plot has been computed in the case (r op t = rmax/f , ro = 2r max /7r, 9 = tt/2). The cartesian 
coordinate [x, z) are such that the axis of revolution of the galaxy corresponds to x — 0. 

r o{@o) such that: <f>o is always the same, and the solutions <J>(r, 0) are exactly identical. We checked this successfully 
on various examples. 

The rotation curve can be deduced from the gravitational potential: 

d 

v 2 = r — $(r, = 7i72). (48) 
or 

On figure 3, we compare the rotation curve obtained following our method with the one calculated in the approximation 
of Paper I: namely, replacing the thin disk by a spherical one, with a density such that in absence of any halo, the 
rotation curve along the stellar plane would be the same as with the true non-spherical disk. In presence of a halo, 
one can see that the difference becomes important only at large radius. 

IV. THE ROTATION CURVES OF DWARF SPIRALS 

Dwarf spiral galaxies are known to be completely dominated by dark matter at all radii. Usual CDM models fail 
to reproduce the rotation curves of those systems. The purpose of this analysis is to investigate whether a self- 
interacting massive scalar field halo is able to reproduce such rotation curves. Therefore, we will first scrutinize the 
typical dwarf spiral galaxy DDO 154 that has been thoroughly studied - see for example the observations by [27] 
and [28]. Because it is isolated and therefore seems to be protected against any external influence, this dwarf spiral 
features a prototypical example for our study. Its HI gas contribution is well measured and follows the distribution 

p sas {r,0) = ^exp{-O.8|cos0| — } exp (-5 1 sin0| — ) . (49) 

I r °pt J I r opt J 
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FIG. 3. Rotation curve due to a disk-shaped density (in red) compared to that obtained in a spherical approximation (in 
green) as in Paper I. These two examples have been computed in the case (r opt = r max /7T , r = 2r opt , 8 = 7r/2). The difference 
is seen to be fairly small. 

Its optical radius r opt is equal to 1.4 kpc. The contribution of its stars is visible and therefore well-known, with a 
density distribution given by relation (45). Both stars and gas account for a small fraction of the observed circular 
velocity. 

In order to compare the various dark matter models with the observations of DDO 154, we have performed a x 2 test 
on the 13 data points from ref. [28]. To commence, we have considered models in which the approximate real density 
of stars and gas has been assumed with /?g as ss 0.15pg tars and u s tars(?*opt) ~ 0.3 - 0.4 Utot(r pt) _ where v s tars is the 
stellar contribution to the rotation velocity «tot- The precise value of the ratio u s tars (?^pt )/ftot (^opt ) is unknown and 
has been adjusted here in order to provide the best fit. On top of stars and gas, a dark matter component is added 
with a density profile that depends on the model at stake. The Moore's model [1] is featured in the panel (a) of fig. 4 
and corresponds to the spherical symmetric density 



where r s is a scale radius parameter that is also adjusted in the fit. Recent CDM N-body simulations point towards 
such a profile. In the case of DDO 154, the best fit corresponds to (ropt/rs) 1 ' 5 ~ 0.07pg tars and to very large 
values of the scale radius r s . The x 2 value is found to be approximatively equal to 600 for 10 degrees of freedom. As 
previously mentioned, Moore's model where the density diverges like r -15 in the central region fails to account for 
the dark matter distribution inside DDO 154. Then, we have tested a NFW spherical density profile [29] where 



Such a distribution peaks at the center and has also been found to naturally arise in N-body numerical simulations 
of neutralino dark matter. We find that the best % 2 lies around 200 when r s w 9 r opt and pJ, FW sa 6 Pg tars - see panel 
(b) of fig. 4. We have also considered an isothermal spherical halo [30] with 




(50) 




(51) 




(52) 
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FIG. 4. From left to right and top to bottom, the best fit of the DDO 154 rotation curve (in red) is featured as well as the 
observed stellar (in green) and gas (in purple) densities. A dark matter component has been added and various profiles have 
been assumed : Moore's density (a), NFW profile (b), isothermal halo (c) and Burkert's phenomenological distribution (d). 



This density has been introduced in order to account for fiat rotation curves. In the case of DDO 154 where the 
circular speed starts to decrease beyond 4.5 kpc, the best fit is obtained for r s « 1.2r opt and pf so « 0.15 p^ taTS . The 
corresponding \ 2 is 110w f ar more better with a value ~ 55 - see panel (c) of fig. 4. Finally, a Burkert spherical 
distribution [31] 

Pb(t) = Pb I j— ° pt 2 2 , > , (53) 

(r + r s ) (r l + rj) J 

has been considered. This density law has a core radius of size r s - just like the isothermal halo - and converges 
at large distances towards a Moore or a NFW profile. It has been introduced as a phenomenological explanation of 
the rotation curves of dwarf galaxies [32]. The best parameters are then r s sa 1.9r opt and p B rs 6/0g tars , leading to a 
best x 2 ~ 45 - see panel (d) of fig. 4. At this stage, we reach the conclusion that neutralino dark matter - should it 
collapse according to the N-body numerical simulations a la Moore or NFW - is too much peaked at the center of 
DDO 154 and does not account for the rotation curve in that region. The fact that these species fail to reproduce 
the inner dynamics of a system known to be saturated by dark matter is definitely a problem. The isothermal and 
Burkert halos provide a better agreement with the data but are not consistent with the decrease observed beyond 4.5 
kpc. 

We have then investigated a slightly different idea. Following Pfenniger and Combes [33], the dark matter inside 
galaxies would consist of pure molecular hydrogen Hi - so cold that it would have gone undetected so far. The 
formation of stars in the inner parts and its concomitant UV light production would have turned part of the H2 
into detectable HI. The distribution of this hidden H2 component could be derived in the case of DDO 154 from its 
observed rotation curve. We will nevertheless adopt the opposite point of view since our aim is to derive - and not to 
start from the circular speed. We have therefore artificially rescaled the observed gas density (49) by a homogeneous 
overall factor. The best fit featured in fig. 5 corresponds to p c gas = Pg tars and leads to a best x 2 °f ~ 500 which is not 
particularly exciting. In the case of the models of fig. 4, the addition of such a cold gas component does not improve 
the goodness of our fits. 
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FIG. 5. The DDO 154 rotation curve (in red) is fitted with the observed stellar density (in green) while the gas distribution 
(in purple) has been artificially enhanced with respect to the observed HI by a rescaling factor. 
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FIG. 6. The best fit of the DDO 154 rotation curve (in red) is presented with the observed stellar density (in green) and gas 
distribution (in purple). In the left panel, a self-interacting bosonic halo is assumed with m 4 /A ~ 75 eV 4 together with the 
observed gas density profile. In the right panel, the gas component has been rescaled in order to improve the goodness of fit 
and a value of m 4 /A ^ 50 eV 4 is derived. 

Finally, we have assumed the presence of a self-interacting bosonic halo and applied the recursion method discussed 
in section III. The left plot of fig. 6 corresponds to stellar and gas populations as observed while a value of to 4 /A « 75 
eV 4 provides a best \ 2 of 16. The agreement with the measured rotation curve is quite good. Notice that the bosonic 
halo dominates completely the inner dynamics beyond ~ 0.5 kpc. More impressive is the right plot of fig. 6 where 
the gas distribution has now been rescaled in order to improve the goodness of fit. A best \ 2 °f ~ 7 is reached for 
Pgas ~ 0-35 Pstars an d a value of to 4 / A ~ 50 eV 4 . 

Beside the prototypical example of DDO 154, we have analyzed a set [32] of small and medium size spiral systems for 
which measurements of the rotation curve are of high quality. These galaxies have been selected on the requirement 
that they have no bulge, very little HI - if any - and a dominant stellar disk that accounts for the dynamics in the 
central region. They are also dominated by dark matter as is clear from fig. 7. A self-interacting bosonic halo has 
been assumed with m 4 /A « 50 eV 4 . Because of the presence of wiggles in the rotation curves - presumably related 
to spiral arms inside the disks - the best \ 2 value becomes meaningless. The qualitative agreement is nevertheless 
correct except in the case of 545-G5 where the optical radius is r opt = 7.7 kpc. Because the mass to and the coupling A 
define a unique scale of ~ 2.3 kpc - see relation (25) - the Bose condensate does not extend enough to account for the 
dark matter inside large systems. A single self-interacting bosonic halo fails to reproduce at the same time the dark 
matter inside light and massive spirals. A possible solution lies in the existence of several small bosonic condensates 
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FIG. 7. From left to right and top to bottom, the best fit (in red) of the rotation curves of various spiral galaxies with 
increasing optical radii. The green lines stand for the stellar contributions. The various panels respectively correspond to (a) 
N7339 (r pt = 4.8 kpc) - (b) M-3-1042 (r opt = 4.8 kpc) - (c) N755 (r opt = 4.8 kpc) - (d) 116-G12 (r opt = 5.4 kpc) - (e) 
563-G14 (r op t = 6.4 kpc) - (f) 545-G5 (r opt = 7.7 kpc). 

or clumps inside the halos of large galaxies whereas a single condensate would account for the dark matter of dwarf 
spirals such as DDO 154. We have concentrated on dwarf spiral galaxies for which ncutralinos seem to be actually in 
trouble. The case of several bosonic clumps is beyond the scope of this work and will be investigated elsewhere. 

V. THE SOLAR SYSTEM 

As long as we were interested in the inner dynamics of galactic systems, the baryonic density p\, in equations (23) 
and (33) was implicitly averaged over distances of order a few pc and behaved smoothly. If the stellar population is 
now made of point-like particles with mass Ms, the gravitational potential $ varies according to 

4 

A$ + 8nG — ($-$o) W($ -$) = 4 7rG^Af l( 5 3 (f-r!) . (54) 
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Outside the Bose condensate, the usual Poisson equation is recovered so that the gravitational attraction of a star - 
say the Sun - is not modified with respect to the conventional situation. Slightly different is the case where the Sun 
lies inside the region where the field <\> extends. Intuitively, the scalar field is expected to be attracted by the solar 
gravity and to concentrate around the Sun whose gravity should consequently be strengthened. In order to investigate 
that effect, we first notice that the potential difference $ — $o is a linear function of the sources within the Bose 
condensate. The contribution $0 of the sun to the potential difference $ — $0 satisfies the modified Poisson equation 

4 

in 

A$ + 8ttG — $0 = 4ttG p Q , (55) 
A 

with the condition that it must vanish on the boundaries of the Bose condensate. If the solar system is well embedded 
inside the latter - below a depth well in excess of a few AU - the surface of the condensate is so far that we may just 
require that $ Q vanishes at infinity. 

Assuming in addition that the solar density p@ has an isotropic distribution, the solution of eq. (55) readily obtains 
in terms of the spherical Bessel functions jo(z) — sinz/z and no(z) = — cosz/ z as explained in Appendix B 

$©(r) = — A-kG — < cosz / pq (ii) 11 sinii du + sinz / p Q (u) u cosu du > . (56) 
r I Jo Jz J 

The dimensionlcss radial coordinate z is defined as the ratio r/L where the typical scale L has already been defined 
in section II. Relations (24) and (25) imply that L exceeds the solar radius i? by some ten to eleven orders of 
magnitude. The gravitational potential which the Sun generates with the help of the scalar field <f> simplifies into 

C M 

d> Q (r) = -— ® cos (r/L) . (57) 

Because of our assumption as regards the boundary condition - which we placed at infinity - this relation may be 
safely used only for distances r <C L. Inside the solar system, this leads to the potential 



and to the gravitational field 



2 V 



G M Q G M q 

ffW = - — " -2£T ■ (59) 



Should the solar system be embedded inside the Bose condensate of the field under scrutiny in this article, the various 
planets and satellites that orbit around the Sun should undergo the additional constant radial attraction 

Sg = 2.85 x 10- 18 cms" 2 {^} {j^Y ■ ( 60 ) 

This acceleration is so weak that it should not alter the motion of the planets around the Sun. The relative increase 
of the solar gravitational attraction is actually 

6g 17 r r l 2 f 1CT 2 ] r m 1 4 

f = L93xl ° {iauI — kv) ■ (61) 



Detailed analysis of radio metric data from Pionner 10 and 11 indicate the existence of an apparent anomalous 
acceleration acting on these spacecrafts [34]. Quite exciting is the observation that this anomalous acceleration Sgp 
is constant and directed towards the Sun. Both features are actually expected in the presence of a self-interacting 
scalar field. However, the magnitude of the observed anomalous acceleration 5gp ~ 8.5 x 10~ 8 cm s~ 2 is ten orders 
of magnitude larger than what is needed to explain the rotation curve of DDO 154. Should the Pioneer acceleration 
be the consequence of a scalar field enhanced solar gravity, it would indicate an exceedingly large value for m/A 1 / 4 of 
order 1.3 keV and a typical condensate size L ~ 9 x 10~ 3 pc. We are therefore lead to the conclusion that we cannot 
explain with the same value of m/A 1 / 4 the Pioneer anomalous acceleration and the rotation curves of dwarf spirals. 
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VI. COSMOLOGICAL BEHAVIOR 



We will now consider briefly the possible cosmological behavior of our scalar field. In Paper II, we studied the 
cosmological evolution of a homogeneous complex scalar field with a quadratic and/or quartic potential. Here, 
we want to update this analysis for the values of m 4 /A found in section IV. Generally speaking, focusing on the 
homogeneous quantities is the first step in any comprehensive study of a given cosmological scenario. In our case, 
we need to know whether the evolution of the field background violates any cosmological bound before studying the 
possible growth of spatial fluctuations - hoping that they will cluster and form galactic halos after the time of equality 
between radiation density and field density. 

We refer the reader to Paper II for a detailed resolution of the Klein-Gordon and Friedman equations in a Universe 
containing ordinary radiation, baryons, a homogeneous complex field and a cosmological constant relevant only today. 
It is straightforward to show that when the potential is dominated by the quartic term, the energy density of the 
field smoothly decays as a -4 : so, in the early Universe, the scalar field behaves as "dark radiation". Later, when the 
quadratic term takes over, i.e., when 



, , 771 

V{4>) ~ 2—, 



(62) 



the field starts to decay as a -3 , like dark matter: so, it can be responsible for a "matter-like" dominated stage. During 
the whole cosmological evolution, the kinetic energy of the field is of the same order of magnitude as its potential 
energy. So, the ratio m 4 /A immediately gives a rough estimate of the total energy density of the field at the time of 
its transition, denoted later as p^ r . If this density is ~ 1 cV 4 , we immediately notice that it is of the same order of 
magnitude as the density at radiation-matter equality - remember that p cq ~ 0.55 cV 4 for the concordance ACDM 
model. So, in the early Universe, the density of our "dark radiation" (the scalar field) had to be comparable to 
that of "true radiation" (photons and neutrinos). This brings some considerable tension with the bound on the total 
radiation density that can be derived from BBN. 

This cosmological toy-model and the problems associated with it were first discussed by Peebles [17] in the case of a 
real scalar field, with essentially the same motivations as in the present work. As a possible way out, Peebles proposed 
a small modification of the scalar potential, 



v(4>) 



m 2 (j) 2 



Xcj) q 



(63) 



where q would be non-integer and slightly smaller than 4. Indeed, by lowering the index q, one can decrease the 
fraction of dark radiation in the early Universe, and in particular at BBN. We will not follow this direction. Indeed, 
the analysis of section IV revealed a preferred value of m 4 /A around 50-75 eV 4 . This is significantly larger than the 
observed value of p cq and than the rough estimate of Paper I where we considered m /A ~ 1 eV . Our purpose in 
the rest of this section will be to check whether this new value is compatible with the BBN bound. 

In order to obtain a precise relation between the parameter m 4 /A and the effective number of neutrinos at BBN - 
which is a convenient way to parameterize a cosmological density which behaves like some extra relativistic degrees of 
freedom we need to study numerically the detailed behavior of the field in a vicinity of the transition between the 
radiation-like and the matter-like regimes. For each value of m 4 /A, it is possible to follow p^, and to extrapolate the 



branches in p$ cx a and in 



oc a 



We define p^I as the energy density given by intersecting the two asymptotes. 



The knowledge of this single number is sufficient in order to relate exactly the constant value of p^a 4 in the early 
Universe to the constant value of p^a 3 measured today. A simple numerical simulation gives 



^transition 



2.4m 4 /A 



(64) 



independently of any other field or cosmological parameters. The simulation also provides a very good analytic 
approximation of the field density at any time - imposing that today, when a = ao, the field density is given by the 
fraction of the critical density usually attributed to Cold Dark Matter, p^ = f2 c d m /9c : 



— O n° 

— i£j cdmP c 



(a /a) 6 + n cdmP ° c 



A 



2.4 m 4 



2/3 



1 1/2 



(a /a) 6 



(65) 



The field density before the transition can be read directly from the previous equation, 
parameterized in terms of an effective neutrino number, defined as usual through 



It can be conveniently 
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AN eS = ^ , (66) 



where p v is the standard density of a single relativistic neutrino species. The final result is 



*=liffifj = 7 -Hi^rJ {—) ■ (67) 

Sticking to m 4 /A = 50 (eV) 4 , and using the currently preferred values h = 0.68 and fi c dm = 0.3, one finds AN c g = 2, 
which is above the usual BBN bound |AiV c ff| < 1 [35]. However, it is still possible to find some values of (AiV e g, h, 
f^cdm) satisfying the above relation, and allowed at the 1-cx level by current CMB experiments and BBN predictions 
[36] - for instance, (1.0, 0.63,0.20). In any case, in a very near future, the new CMB observations will set some 
stringent limits on these three parameters: it will then be easy to state about the validity of our alternative to the 
usual cosmological scenario. 



VII. CONCLUSIONS 



We have shown that a self-coupled charged scalar field provides an excellent fit to the rotation curve of the dwarf 
spiral DDO 154. That galaxy is the prototypical example of a system known to be completely dominated by dark 
matter. The effect of the quartic coupling A results into an effective modified gravitation inside the Bose condensate 
where the Poisson equation becomes strongly non-linear. The problem - complicated by the non-sphericity of the 
baryon distribution - has been solved exactly as explained in section III. The agreement with the observations of the 
circular speed of DDO 154 is impressive. Notice that neutralino dark matter does not pass this test because of the 
central cusp that it would develop. We conclude that the charged scalar field considered in this analysis provides an 
exciting alternative to the galactic dark matter - at least inside dwarf systems. A typical value of m 4 /A ^50-75 
cV 4 obtains. 

The scalar field behaves cosmologically as a dark radiation component as long as the quartic contribution of the 
potential V (</>) dominates over its quadratic counterpart. The situation gets reversed when the field energy density 
is ~ 2.4 m 4 /A and a matter-like behavior subsequently ensues. The larger the crucial parameter m 4 /A, the sooner 
the transition between dark radiation and dark matter-like behaviors and consequently the smaller the contribution 
of the scalar field to the overall radiation density at early times - for a fixed scalar field mass density today. A large 
value of m 4 /A translates into a small number of effective neutrino families during BBN and we have shown that our 
model marginally satisfies the requirement that AiV e ff should not exceed 1. 

Actually the model is strongly constrained on the one hand side by the size L of the Bose condensates - and therefore 
of the corresponding galactic halos - and on the other hand side by the contribution AiVeff to the radiation density 
at BBN. Both L and AiVeff decrease with m 4 /A and a value for the latter of ~ 50 - 75 eV 4 - which provides excellent 
agreement with DDO 154 - is marginally consistent with BBN. Large halos cannot consequently been pictured in 
terms of a single Bose condensate and the simple scheme presented here has to be modified. A possible solution - 
yet not very natural - is to replace the quartic field self-interaction by a <p q term as suggested by [17]. This would 
alleviate the BBN constraint. 

Another option worth being explored is to imagine that massive and extended halos are formed of several bosonic 
clumps. The coherent configuration that has been investigated here may be understood as the ground state of some 
gigantic bosonic atom. It is therefore conceivable that the scalar field may also form several such condensates that 
would be organized inside a huge bosonic molecule with a spatial extension much in excess of L. The electron cloud 
around the proton docs not extend further than ~ 10 -10 m inside the hydrogen atom and yet electrons are delocalizcd 
over meter size distances inside metals. 

If so, the dark matter would be made of small bosonic clumps. Should the solar system lie within such a system, 
the motion of its planets would provide in that case a lower bound on L since the smaller is the latter, the stronger 
is the effective modification to Newton's law of gravitation. We have actually shown in section V that the scalar 
field concentrates in the solar potential well and strengthens it to generate an additional gravitational attraction that 
is radial and constant. As a matter of fact, the radio data from the Pioneer probes are consistent with such an 
anomalous acceleration that seems to be constant and directed towards the Sun. Assuming that it results from the 
self interacting scalar field which we have investigated in this work, the observed magnitude 6gp ~ 8.5 x 10 -8 cm s~ 2 
would imply a value for m/A 1 / 4 of order 1.3 keV and a typical condensate size L ~ 0.01 pc. The merging of many of 
these small bosonic clumps into a larger structure like a galactic halo is an open question. 
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APPENDIX A 

The purpose of this section is to derive the gravitational potential generated by a pure scalar field condensate. This 
calculation is complementary to the one at the end of section II, based on the polytropic formalism. In absence of a 
baryonic density pb, eq. (23) reads 
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A ($ - $ ) + 87rG— ($ - $ ) H ($o - $) = . (Al) 
A 

Since there is no source, let us suppose that the gravitational potential has a spherical symmetry. One can do a 
change of variable r = Z \J 8v q,^ so that the equation becomes simply 

d 2 

— {z ($ - $ )} + *($ - *o) W($o-*) = . (A2) 

The only solution of this equation which is continuous and derivable everywhere, and goes to zero at infinity is 

sin z 

$o + ($(0) - $0) if z < 7T (see Paper II) 

7r$ (A3) 
otherwise. 



So. the maximum extension of the scalar field halo is r max = n\l - — — — r . The density of the bosonic halo is 

" 87tGto 4 



A 

and the total mass is 

3/2 



2ra 4 

(*-$o)W(«-7r) (A4) 



\ \ j / ■* r+00 _ / \ 

A \ / , ,2 , , 7T / A 



On the other hand, the Gauss theorem applied to the sphere of radius r = r max gives 



GM GM 8nGm 4 
*o = = V - T - A6 

'max 7T V A 



so that $(0) = 2$o- The solution (A3) can be rewritten as 



$o(l + — — ) if z < I" 

(A7) 

otherwise. 



APPENDIX B 



The purpose of this Appendix is to provide a prescription for the definition of - the starting function in the 
recursive method. We tested this prescription on various examples, and found that is always a fairly good 
approximation of the exact solution, allowing for quick convergence. 

The idea is to enforce the boundary surface on which the field density (m 4 /A)($ — <I>o) H(^o~ vanishes to be 
a perfect sphere. Of course, this has to be wrong when the baryonic density is non-spherical. So, if we impose a 
constant boundary radius, we need to relax the fact that the value of $ should be constant all over the boundary. In 
other terms, if r is the radial coordinate, we replace 7i($o — *&) by H(ra — r), where ro is an arbitrary boundary radius. 

If we define a dimcnsionlcss radial coordinate z = r^J^^-, eq. (23) becomes 

A<f> = S(z,6) + ($ -$)H(z -z). (Bl) 
where S = AitGp. After a Legendre transformation, we obtain the following set of differential equations: 

~ (z 2 ^) - ^ (*i - M«,o)*(*b -z) = S,(z) (B2) 



18 



where 5i : o is the Kronecker symbol. Each of these equations can be solved separately on the two intervals < z < zo 
and zq < z < oo, using Green functions. In terms of the spherical Bcsscl functions ji(z) and ni(z), the solution of 
eq. (Bl) for z < zq is 

§i(z) = $ Q S l0 +n l (z) / u 2 j l (u)Si(u)du+ ji(z) u 2 m(u)Si{u)du + C}ji{z) (B3) 



where C\ is a free constant. In the same way, for z > zq, the solution is 

* l{z ) = --L- Z -W [ Z u ^ Sl {u)du--±-z l l^u^S^du-Cf-^-z-^ (B4) 



21+1 J ZQ ,w 2Z + 1 7 Z lw 1 2/ + 1 

where C 2 is another free constant. We impose that each Legendre coefficient of the gravitational potential is derivable 
and continuous on the boundary z — zq: 



This defines a unique value for each constant of integration: 

c x = (J?-J t 1 -<Mi,<))(i + l) + zo(./, 2 - J, 1 ) 

' (/ + i)ii(^o) + ^oiK^o) 

2 _ -z<+ 2 (2; + l)j,(z„)( J? - J I) + z l + 2 {2l + l)j>izo)(I? - If - $o S lfi ) 



(B5) 



(B6) 



where 



C 

1 (I + llnizo) + zqjKzo) 



Ij = ni(z ) I u 2 ji(u)Si{u)du. 



1 



21 + 1 



rz 

J i = n'i(zo) I u 2 ji(u)S[(u)du, 



(B7) 



/ 



J? = / u^S^du. 



One can then reconstruct the gravitational potential from 



+ OC 



&°)(z,0) = ^P;(cos6i)$;(z). (B8) 



1=0 



So far, the approximate solution $' ' constructed in this way depends on two arbitrary numbers: first, Zq, and 
second, $o, which appears explicitly in the definition of Cq and Cq. However, zq and <&o have to be related in 
some way. Indeed, if $(°) was an exact solution, <j> would be equal by definition to &°\zo,6) for any 9. In our 
approximation scheme, <E>(°' (zoi is not independent of (9, but we can choose a particular direction 9q, and impose 
that <&o = ^ ^{zo,9a). Inserting this identity in eq. (B8), and using cq. (B3), one obtains the relation 

$o = <f> {0) (z Q ,9o) = t^—Y' Pi(cos9 ) (ni(z ) f ° u 2 ji(u)Si(u)du + K Ul (z Q )\ (B9) 

where 

(I 2 -ll)(l + l) + z (J 2 -Jl) 
1 (l + l)j l (z ) + z j' l (z ) ■ 1 UJ 

In summary, the first step of our recursive method is performed in the following order: 
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1. we choose a value zq (or tq) and an arbitrary direction 6*o (that will be kept for all the following iterations). 

2. we solve eq. (B9) in order to find $o- 

3. we compute $^(z,9) using eqs. (B3), (B4), (B8). 

The next iterations are performed in the much simpler way described in section III. 
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